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Abstract -The mutual interaction between Cooper pairs is proposed as a mechanism for the 
superconducting state. Above T c , pre-existing but fluctuating Cooper pairs give rise to the un¬ 
conventional pseudogap (PG) state, well-characterized by experiment. At the critical temperature, 
the pair-pair interaction induces a Bose-like condensation of these preformed pairs leading to the 
superconducting (SC) state. Below T c , both the condensation energy and the pair-pair interaction 
0 are proportional to the condensate density N oc (T ), whereas the usual Fermi-level spectral gap 
A p is independent of temperature. The new order parameter d(T), can be followed as a function 
of temperature, carrier concentration and disorder - i.e. the phase diagrams. The complexity 
of the cuprates, revealed by the large number of parameters, is a consequence of the coupling of 
quasiparticles to Cooper-pair excitations. The latter interpretation is strongly supported by the 
observed quasiparticle spectral function. 


I. Introduction. — As Occam’s razor would suggest, 
second order phase transitions often depend on few pa¬ 
rameters 1 . Such is the case for the familiar mag¬ 
netic, spin glass, charge-density wave, structural transi¬ 
tions, etc., and the more exotic Kosterlitz-Thouless case 
for two-dimensional systems 2 . The superconducting 
phase of ‘classical’ materials follows this trend wherein 
a weak attractive electron-electron interaction is respon¬ 
sible for the transition. Essentially a single energy scale, 
the SC gap parameter Ao at zero temperature, is rele¬ 
vant [3]. In the conventional theory of Bardeen, Cooper, 
Schrieffer (BCS) [ 3 ] pair-breaking quasiparticle excitations 
restore the normal-metal state at the transition, such that 
A 0 = 1.76 kBT c . Moreover, A 0 fixes the scale of the crit¬ 
ical currents, for example in a Josephson junction, the 
upper critical field and the coherence length, such as the 
vortex core radius. 

In highly-correlated electron systems such as cuprates, 
many physical properties of their SC state are not yet fully 
understood. In addition to many theoretical ideas (4j|9 
the experiments reveal an inherent complexity: a complex 


quasiparticle self-energy 
lective mode 


10-12 


13181, anomalous specific heat 19 


coupling to a spin col- 
mult i- 


ple energy scales [20-23 , spatial inhomogeneities 24- 29 
checkerboard oscillations 30,31 , pseudogap phenomen; 
above T c [7|[l9[[32f 35], etc. The SC state of cuprates i 
indeed far from the desired lex parsimoniae. 


The clear jump in magnitude of the physical parameters 
is also striking. In addition to T c , the very small coher¬ 
ence length, large penetration depth, large Fermi-level gap 
value, small lower critical-field, etc., come to mind. In the 
case of cuprates, the precise link of the measured Fermi- 
level gap to the condensate is still controversial. Consider 
BiSrCaCuO (2212) as a typical example. The energy gap 
is A 32 meV near optimal doping (p) which, if directly 
tied to the critical temperature using the BCS ratio, would 
give the erroneous value of T c ~ 220 K. Moreover, on the 
underdoped side of the phase diagram the energy gap in¬ 
creases while T c decreases - a clear paradox. Thus, A is 
not the scale of T c . 


The quasiparticle (QP) spectral function, expressing the 
material’s microscopic interactions, also contradicts the 
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conventional picture 10 -12 18,36. Measured by tunnel¬ 


ing or ARPES as a function of rising temperature, the 
cuprate SC gap does not close at T c as would be expected, 
but rather a pseudogap remains at the Fermi level. To 
within thermal broadening, the pseudogap has about the 
same value as the SC gap, with no coherence peaks, and it 
finally vanishes at the higher temperature ^ T* 37 . An 


analogous PG exists in the vortex core 38 -40 or due to 
disorder 25,28 , where SC coherence is also lost. Thus in 


the case of cuprates, the direct sign of the condensate is 
not the gap per se , but instead the finer structure of the 
QP spectral function (peak-dip-hump). 

Very different interpretations have been proposed: the 
two-energy scale gap function, implying a non-retarded 

and the coupling to a spin 


pairing interaction 
collective mode 



Although phenomenologi¬ 


cal, the static model describes very well the particular sig¬ 
natures of the spectral function, the SC density of states 
(DOS), and the transition to the pseudogap state 22 


Still, the origin of these signatures remains to be clarified. 

Many studies with scanning tunneling spectroscopy 
(STS) (24}§9||40 , showed real-space inhomogeneities even 
for pristine samples at low temperature. According to con¬ 
ventional wisdom, without an applied magnetic field or 
current, the SC order parameter should be homogeneous. 
The phenomenon has led to in-depth theoretical studies 
of the mechanisms underlying long-range SC order with 
particular emphasis on the effects of disorder and perco¬ 
lation 47-50 . By studying a variety of spectral data, we 
find that they have a distinctive shape on the mesoscopic 
scale; the underlying parameters then become meaningful. 

To summarize, a solution of the high-T c problem must 
account for: 

( i ) multiple energy scales and large spectral gap A p , 

(ii) large T c , but paradoxically: VbT c <C A p , 

(Hi) the dome-shape of the T — p phase diagram, 

(iv) the SC to pseudogap transition, 

( v ) the unconventional QP dispersion (peak-dip-hump). 

In this work a microscopic theory is presented that ad¬ 
dresses these requirements. We assume the existence of 
pre-existing fluctuating pairs above T c in the incoherent 
pseudogap state [6}f9, 21,34|f5l]. Since the carrier-density 
and coherence length are small, a system of interacting 
bosons is relevant; a new pair-pair interaction, /%j, leads 
to a Bose-like condensation to the coherent SC state at 
the critical temperature (section III). 

Simple expressions for the gap equation and the 
condensation energy are given and the nature of the SC 
to PG transition is described. A new order parameter 
/3(T) naturally emerges, whose properties are examined 
as a function of temperature, carrier concentration and 
disorder - i.e. the phase diagrams. It extends the 
notion of ‘pair field’ that we reported in the coarse-grain 
approach 42 , and gives a direct link between the SC 


coherence and the parameters of the QP spectral function. 


While the Bose-like condensation of preformed pairs is 
relatively straightforward, the interpretation of the exper¬ 
imental QP spectra suggests a much larger number of pa¬ 
rameters. In this work we argue that the further complex¬ 
ity is due to the coupling of quasiparticles to Cooper-pair 
excitations. This proposition is strongly supported by the 
variety of experimental spectra, analysed in sections IV 
and V, and other widely accepted properties of high-T c . 


II. Hamiltonian for interacting pairs. — In the ab¬ 
sence of pre-formed pairs, the so-called ‘pairing’ Hamilto- 
takes the form: 


man 
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H = H 0 + H 


pair 


H 


(1) 


where, in standard notation, 


H 0 = Y j 

k,cr 


£k Ck,cr 


t 


Qc,cr 


( 2 ) 


Hpair E (Afe&fet + Afc'&fc-Afe <6* >) (3) 

k 

with: b k = c k, t an d H is bilinear in the fluctu¬ 

ation terms: bk — < bk >. Now, in the problem we are 
addressing, we consider pre-existing pairs that we label by 
i. We assume that these pairs are still described by the 
BCS model, so that eq. © takes the form: 

iw = -EE( A * 6 * t + Ap&i) (4) 

i k 

after dropping the constant term. We assume that the 
bilinear fluctuating terms in H that are diagonal in i only 
renormalize the pair amplitude. On the contrary, we keep 
the non-diagonal terms such that a new mutual interaction 
between pairs emerges. We find H = H int can be written : 

Hint = \'Z'Ztii'K.bi i + h.c. (5) 

i,j k,k' 


where f3 k ° k/ is the microscopic coupling. This term is ne¬ 
glected in the conventional approach. 

In the absence of the pair-pair interaction Hi nt = 0, 
then Hpc = Ho + H pair describes a pseudogap state of un¬ 
condensed pairs having a binding energy compatible with 
~ T* ^ T c . As expressed in eq. 0. they have a fluctu¬ 
ating amplitude distributed according to: 


Po( Al) 


°o 


(A*. - A 0j fe) 2 + ctq 


(6) 


where A 0 ,/c and ao are the average pairing energy and 
width, respectively. These two parameters will remain 
central to our model throughout this work. 

In the non superconducting state, the excitation spec¬ 
trum is: E l k = + | A^| 2 with the corresponding spec¬ 

tral function: 


A(e k ,E) = -- Im 

7r 


f 

Jo 


dAi Po(Ai) 
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E - El + i r 


(7) 


P-2 













































Pair-pair interactions 


where T is the Dynes QP broadening parameter 53 . As 


shown in 42 , the above spectral function and correspond¬ 
ing density of states (DOS) leads to a smeared gap of width 
Ao at the Fermi level characterized by vanishing coherence 
peaks (similar to Fig.[4j spectrum 5). 

The interaction between Cooper pairs, P k 3 k /, is respon¬ 
sible for the superconducting state. As in our previous 
work 42 , this interaction is of the form: 


^=9k9k'P 0 (Ai)p 0 (Ai,) 


(8) 


where the factor g k preserves the d-wave symmetry, and 
Po(Afc) is the pair distribution. With this potential the 
equalizing of the A^. in the final SC state is favored while 
retaining the memory of the initial state. It has the useful 
property of being separable, which allows for the decou¬ 
pling of the equations. 

In the mean-field approximation, the operator b\, is 
replaced by its quantum average: < b k , >= A J k ,/(2E k ,), 
so that the effective interaction reduces to: 



Temperature [Kelvin] 

Fig. 1 : Phase diagram in the interacting-pair model. 

with the final gap equation: 

A fe = Ao, fe -2^P 0 (A fe ) (15) 


H int = J2 2 ^Po(^l)b^ + h. 


i,k 


where, /?fe = 9k E 9k ' P o(A {,) 


j,k' 


At 

8 El 


(9) 

( 10 ) 


which defines p k in terms of all other pairs j 7 ^ i in the 
system. 

We assume that in the PG state, 3h. is negligible due 
to the fluctuations of A^, see eq. (10). On the contrary, 
upon condensation, A J k —>• A&, one has: 

Pk = Ap x g k ^ g k > P 0 (A k t) 


A k > 

$E k , 


( 11 ) 


which is much larger due to the factor N p , the number of 
pairs. Moreover, eq. is a new self-consistent equation 
which depends on the exact A k and the corresponding QP 
excitations : 

E k = \/ e l + Af 

To obtain the SC gap equation, we first assume that all 
pairs join the condensate (a condition we relax in Sec. Ill): 




b k P Then H pa i r takes the form: 




pair 


■53(Ao,fc6fc t + A 0 , fc *6 fc ) (12) 


where A 0 ,/c is average value in the initial state ((3 k = 0 ). 

The interaction term Hi nt is evaluated directly in the fi¬ 
nal state (A k —)> A*;), in the spirit of the Brillouin-Wigner 
perturbation approach: 


Hint = 2 a(Afe) b k ] + h.c. (13) 

k 


Combining (12) and ( [l3| ), the total Hamiltonian is then: 


H = H 0 ~Yi (Afc & 4 +A fc * 6 fc ) (14) 

k 


Written in this form, the equivalence of the separable pair- 
pair interaction | 8 | to an additional ‘field’ is evident. 

In this work, A k is assumed to be d-wave and, at the 
Fermi-level (e k = 0) along the anti-nodal direction, we 
simply note A k = A p : 

A p = A 0 -2pP 0 (A p ) (16) 


We stress that the effect of the interaction (- f3 above) is 
to lower the pair binding energy with respect to Aq. 

The maximum (anti-nodal) gap A p appears on both 
sides of eq. (16) which must be solved self-consistently. 
Defining the change in energy at T = 0 as 


s c — Ao A p 

we obtain the cubic equation in £ c : 

2 

e c = 2/3P 0 (A p ) = 2/3^A_ (i 7 ) 

£ c -T <7 0 

which shows that £ c has a strong dependence on the in¬ 
teraction P as well as the distribution width, ao. For 
optimally-doped BiSrCaCuO(2212), £ c ~ 22 meV, which 
is still large compared to fc^T c , but less than the spectral 
gap A p ~ 38 meV. The three characteristic energies: A p , 
e c , and A 0 = A p + £ c , are illustrated in the phase diagram 

Fig-0 

The interaction P has a dominant effect on the energy 
changes involved in the transition. To illustrate, consider 
the BCS expression for the self-consistent gap: 

A p = 2 E 0 e~ 1/iNoV ^ 

where Eq is the cut-off energy, Nq is the Fermi-level DOS 
and V p is the total pair potential. Using V p = Vo + Vint = 
Vo(l —/?), where P is the relative pair-pair interaction, and 
Ao and £ c are of the form: 

A 0 = A p (1 + P) and e c = pA p (18) 
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Fig. 2: Influence of a pair-pair interaction, (3 

As illustrated in Fig. [2j the spectral gap A p decreases 
sharply as a function of /3, becoming vanishingly small as 
/3 1 . Ao begins equal to A p , for /3 = 0, and then follows 

a similar trend as A p , albeit with Aq > A p . The main 
observation is that £ c is nearly dome-shaped - a general 
property of high-T c and indeed other types of supercon¬ 
ductors. Here, the pair-pair interaction f3 causes the rapid 
decrease in A p and the dome-shape of the condensation 
energy. 

The physical origin of the pair-pair interaction within 
the condensate can now be addressed. We consider that 
it is a non-retarded interaction mediated directly by the 
quasiparticles. Consequently it is long-range, implying 
that f3 is proportional to the number of condensed pairs 
N oc (T ). It can therefore be written: 

0{T)=p o N oc (T) (19) 


The effect of this interaction on the condensation mech¬ 
anism is investigated in the following Sec. Ill and on the 
quasiparticle spectrum in Sec. IV. 


III. Pair condensation. — Consider the phase dia¬ 
gram of Fig. [I] in the context of uncondensed pairs above 
T c . Qualitatively, the total pair energy A p follows the crit¬ 
ical curve shown and which vanishes at ~ T*. In particu¬ 
lar, the SC critical temperature, denoted by C, is situated 
on the plateau of the critical curve where A p is about equal 
to its zero-temperature value. The fact that A p remains 
constant in the SC to PG transition can be directly seen 
as a function of temperature 36 37 , disorder 
in the vortex core 


38 39 


25,28 and 


In the BCS gap equation, the reduction of the gap with 
rising temperature below T c is associated with quasiparti¬ 
cle excitations given by the factor : 


1 - 2 f(E k ,T) = tanh 


Eh 


2k B T 


where f(Ek,T) is the Fermi-Dirac function and E^ the 
QP energy. In the case of cuprates, as the temperature 
approaches T c from below, the smallest value of this func¬ 
tion is tanh(A p /(2 k B T c )) ~ .96 and is ~ 1 otherwise. 


Clearly, pair-breaking linked to thermally excited quasi¬ 
particles is not relevant. 

Assuming 2D free-electrons for a single conducting 
plane, we note that the pair density per unit area is about: 
N p /A ~ mA p /(2irh 2 ). Then, consider the typical number 
of pairs within the disk of area S = i r£ 2 , with ^ £ being a 
pair diameter: 


N p ~ 


mAp £ 2 

2 h 2 


~ 1-2 


Clearly, the combination of small £ and small pair overlap 
strongly supports an interacting boson model. 


Table 1: Summary of main parameters 



A p 

self-consistent SC gap 

SC gap equation 

p 

pair-pair interaction 


Ao 

pair distribution average 


^0 

pair distribution width 

Condensate 


excitation mini-gap 

QP spectrum 

a 

pair broadening 

r 

lifetime broadening 


We thus retain the preformed-pair hypothesis above T c 
and propose that their condensation follows Bose statis¬ 
tics. To proceed, one needs the excited states of the system 
which, in our model, are given by the previous distribu¬ 
tion Pq(AQ. In this Section, the magnitude of the gap 
is relevant to these excited pairs and we drop the explicit 
^-dependence. 

The occupation number, 7V oc (T), vanishes in the non¬ 
superconducting PG state, N oc (T c ) = 0, all pairs being 
in excited states. On the contrary, 7V OC (0) = 1 in the 
SC state at zero-temperature, at which all pairs have the 
unique value: A 2 = A p . Analytically, the model would 
then imply: 

POO 

N oc (T ) = A dA* P 0 (AQ f B ( A* - A P ,T) (20) 

J A p 


f B (A l ,T) is the Bose distribution, with vanishing chem¬ 
ical potential {fi B 0), A 1 — A p is the pair excitation 
energy and A is for normalization. It is evident that this 
integral diverges since at the lower limit where A 1 —>■ A p 
the distribution Pq(A 2 ) remains finite while f B (A l ,T) di¬ 
verges. This well-known singularity indicates that a 2D 
condensation of non-interacting bosons cannot occur [54] . 
Furthermore, a Kosterlitz-Thouless transition can also be 
ruled out at this temperature scale. 

There must therefore be a strong attractive interaction 
between a given boson and the condensate. Since the Bose 
distribution in eq. (20) decreases sharply with energy, the 
low-lying excitations just above A p are critical. A sim¬ 
ple model for the energy needed to remove a pair from 
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Fig. 3: Temperature dependance of the Bose condensate pair density N oc (T) using eq. (21). 

Upper right: Fit to QP DOS peak height from [57] (reported T c = 83 K) 

Lower right: Josephson I c Rn from SIS break-junction [55] (reported T c = 77 K) 

Left panel: N oc (T) corresponding to 5 samples from overdoped (1) to underdoped sides (5). 


the condensate is a mini-gap S separating the condensate 
energy from the excited states. We thus write: 

POO 

N oc (T) =A dA* P 0 (A i ) f B (A* - A p , T) (21) 

to calculate the occupation number. The free-energy at 
the transition can likewise be written 


54 


P{T C ) = A 


I 

J a 


Ap+<5 


dA* PoiA*) giA* - A P ,T C ) (22) 


where g(A\T) = —ks Tln(l — exp(— A 1 /(ksT)). 

The interpretation of the mini-gap S will be discussed 
later. In practise, we use the experimental value of T c to 
determine the normalization A of the pair distribution. 
For an optimally-doped BiSrCaCuO(2212), implementing 
eq. (21) with S = 2.1 meV and T c = 95 K leads to E(T C ) ~ 


0.84 ks T c , which is a satisfactory order of magnitude. 

The above Fig. [3] summarizes just how well the Bose 
condensation model fits the data. First, precise fits to the 
QP tunneling spectra from Renner et al. 37 , as a function 


of T, allow to determine the parameters A 0 and do, as well 
as the QP peak heights (PH). As we showed in our previ¬ 
ous work [22,42 , the attenuation of the DOS peak height, 
analogous to the effect shown in the spectra of Fig. [4] is 
due to the gradual decrease in the interaction /3(T), as 
expected from eq. ©. The condensation energy e c {T) is 
also proportional N oc (T ), eq. <EZf due to its dependance 
on /3(T). 

In the upper right panel, we plot both the data points 
from Renner et al. 37 and the best fit of N oc (T ) using 


eq. (21), which is excellent except for a slight jump near 
T c . Unfortunately, similar STS experiments as a function 


of temperature are rare. In the lower right panel of Fig. [3] 
we plot the critical current I c Rn obtained by Miyakawa 
et al. using a SIS break-junction setup [55]. We empha¬ 
size that the Josephson I c Rn is directly sensitive to the 
SC condensate. Here it decreases monotonically with T 
whereas the gap value A p remains constant in this tem¬ 
perature range, in agreement with our model. In both 
experiments, there is a shift in T c comparing the local 
vanishing of N oc (T ) and the reported bulk value. Still, 
the quantity /3(T) oc N oc (T ) has the required properties 
for an order parameter. 

The overall behavior of N oc (T ) from underdoped to 
overdoped BiSrCaCuO(2212) is illustrated in the main 
panel of Fig. [3] The curves are quite different from the 
conventional BCS order parameter having a short plateau 
towards low-T (spanning ~ 10 K) then descending rapidly 
to intersect the T-axis with a finite slope. Starting from 
the fitted estimate of the mini-gap, 5 = 2 meV, we scaled S 
with T c using 3.85 J = T c . The values of the parameters 
A 0 and cr 0 of the pair distribution were obtained from pre¬ 
cise fits to the experimental spectra on BiSrCaCuO(2212), 
the subject of the following section. 


IV. Quasiparticle spectrum. — The QP spectral 
function and DOS, obtained by ARPES and tunneling re¬ 
spectively, contain valuable information on the SC state 
at the microscopic scale. The SC spectral function is: 

A(e k ,E) = -~ Im 1 (23) 

7 r h — Rk + ^ 1 


where is the QP energy. The BCS-Bogoliobov coher¬ 
ence factors are omitted since, as described by Schrieffer 


56 , they are redundant if A/ C (±e/ C ) is symmetric with re¬ 


spect to the Fermi level. 
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Fig. 4: Left panel: Fits to a series of 5 QP spectra from McElroy et al. 28 using eq. (28) (dots: experiment, solid line: 


theory). The spectra are ordered according to the ‘strength’ of their SC features. Insert: Experimental spectra on the 
occupied side before removing of background slope. 

Right panel: Values of important parameters from the spectra. A 0 , cr 0 , and Edi V are all moving upwards, following A p . 
The pair-pair interaction /3 oc N oc moves oppositely towards lower energy, parallel to the peak heights (PH). 
a varies from 2 to 16 meV from spectrum 1 to 5 (not shown). 


The QP DOS is obtained directly from the spectral dependant gap: 
function using: Ns(E) = 2 k A(e k ,E) where, for a d- 
wave pairing symmetry, A k (6) = A p cos(2 0) leading to the 
typical V-shaped DOS 40, |57||58] . However, as we have 


N k (E k ) = A 0 ,/e -2 /3 k Po(E k ) 


(25) 


shown repeatedly 22 ||41||42j , such a d-wave DOS, with 
weak logarithmic singularities at the coherence peaks and 
no dip-hump structures, fails to match the experimental 
spectra. 

We have shown that the total Hamiltonian 0 with the 
interaction term 0 not only expresses the ground-state, 
with all A l k within the condensate A^ = A&, but also the 
pair excited states, with the corresponding distribution, 
Po(A k ). The latter excited states lower the condensate 


occupation number N oc (T ), eq. (21), and eventually the 
SC coherence is lost at T c . How does this relatively simple 
Bose-like transition affect the QP spectral function ? 

In our model, the pair-pair interaction is directly medi- 


Note that the QP dispersion, while having the same func¬ 
tional form: 

E k = (26) 

is in fact greatly modified due to this energy dependence 
and must be used for calculating the corresponding DOS. 
With surprisingly few parameters, this QP dispersion and 
gap function eqs. ( 25|26 ) accurately matches the wide 
peaks and characteristic dips as measured in the spectra, 
see Fig§ 

The derivation of the DOS is given a full treatment 
in j4l][42 . We recall that, neglecting T-broadening, it 
depends on the following derivative: 


m our model, tne pair-pair interaction is directly medi- jy 0 r°° 

ated by the quasiparticles. An excited pair, i.e. A\ > A&, Ns(E,6) = — J ^ de k S(E k — E) 


can be coupled to a quasiparticle of the same energy and 
conversely, that an excited QP of energy E k is linked 
to such a ‘virtual’ pair. We propose that this energy- 
conserving coupling is manifested directly in the pair-pair 
interaction term of the Hamiltonian, eq. §. Indeed, re¬ 
placing Al by El we have: 


No 

2 7r 


de k 


9E k 


H int = 2p k P 0 (E i k )bi i + h.c. 

i,k 


fe J E k =E 

(27) 

where Ns(E , 0) is the partial DOS in the 0 direction. The 
final DOS therefore involves the energy-derivative of the 
gap function itself - a rather unusual property. A full 
treatment, using the spectral function, leads to: 

2w d9 E — iT — A k (E, 9) aA g^ ,6>) 


(24) 


N S (E) =N 0 Re — 


l 


2 ?r y/(E-i r) 2 -A fc (£,0) 2 


In the final SC state, this leads to a non-retarded energy- 


(28) 

This unique DOS function is used to fit all the spectra 
in this work, Figs j4] and [5| where the spectral shape is 
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determined by the distribution Po(^), i.e. Ao and ao, 
and the value of ft (see Table [I]). 

Above T c there is the additional effect of thermal pair 
excitations above A p . Also, pair excitations arise due to 
disorder, which we take into account by the broadening 
parameter a of the pair potential: A p -» A p — i a in 
eq. (28). Its value is important when N oc (T ) < 1 or at 
low-doping (as shown in Fig. [6|. Otherwise, in all fits, 
r r ^ j 1.5 meV, is small. 

full 


In the detailed experiments by McElroy et al. 28 


dip position E&i p . On the contrary, the pair interaction 
parameter fi moves in the opposite direction , concomitant 
with the lowering of spectral peak heights (PH). 

This observed increase of gap width A p with decreas¬ 
ing p is a general consequence of the pair-pair interaction 
model, eq. (16). Moreover, the discussion on the temper¬ 


STS mappings of BiSrCaCuO(2212) have been performed 
to study the local SC characteristics and, with atomic res¬ 
olution, the effects of the local oxygen dopants. The study 
of three different samples with nominal average gaps (45, 
55 and 65 meV) are of particular interest, ranging from 
near optimal to increasingly under-doped. The DOS vari¬ 
ation due to atomic oxygen dopants is quite significant: 
a resonance near -.96 meV is seen on the occupied side 
of the spectrum and leads to a strong variation of the SC 
spectral shape. 

Two other points in these experiments are significant 
for the present model. First, local spectra from one sam¬ 
ple can be identified with spectra from another sample at 
locations having common characteristics (gap width, peak 
height, dip position). Although the samples show topolog¬ 
ical variations, this suggests that on the mesoscopic scale, 
the SC gap function has consistent properties independent 
of the sample and its parameters are physically meaning¬ 
ful. Secondly, the optimally-doped sample is much more 
homogeneous than the other two - suggesting that a per¬ 
colation effect must be taken into account for very under¬ 
doped samples. 

The spectra in Figj4] are sequenced according to the 
‘strength’ of their superconducting characteristics: in 
spectrum 1., a smaller gap width, with sharper QP peaks 
and in spectrum 5., a larger gap and very attenuated 
peaks. The spectra with strong SC features largely dom¬ 
inate on the optimally-doped sample, whereas the weak 
SC spectra are found in mesoscopic regions in the inho¬ 
mogeneous under-doped sample. The series resembles the 
SC to PG transition at low-temperature as seen by Ren¬ 
ner and Cren [25,37 and in this context we consider the 
compatibility of our model. 

In Fig j4j left panel, we show the fits to the spectra, 
having removed the background slope and symmetrizing 
in each case. In the right panel we plot the values of 
the parameters in the DOS function eq. (28) and their 
variation from spectrum to spectrum. 

The first point is that our gap value A p is 10-20 % 
smaller than the nominal ones cited by McElroy et al. In¬ 
deed, the observed spectral gap is affected by the deriva¬ 
tive of the gap function at — A p , as illustrated in Fig [5] 
lower panel. The measured peak is also affected by the pa¬ 
rameter a, which reaches ^ 16 meV for spectrum 5. This 
damping effect is smaller in the case of sharp QP spectra. 

Secondly, it is remarkable how well several parameters 
scale with the spectral gap, in particular Ao, cr 0 and the 


ature dependance of /3, where we have /3(T) oc 7V oc (T), 
leads to the conclusion that the condensate occupation 
density N oc is decreasing from spectrum 1 to 5. This is 
confirmed by the quantity £ c = Ao — A p , which remains 
relatively robust in this sequence, and <jq which becomes 
large. The likely scenario is the effects of disorder in the 


underdoped regime 47,50 


In summary, the simple QP spectrum, eq. (28), matches 


remarkably well the measured local DOS of an inhomoge¬ 
neous SC state, e.g. spectra 3-5. In a given region of the 
sample, the main effect is thus the change in magnitude 
of the condensate density N oc (T ) on the mesoscopic scale 
and, through the pair-pair interaction /3, the damping ef¬ 
fect on the QP features (peak-height and dip strength). 
Lastly, spectrum 5 shows very weak quasiparticle peaks, 
close to the PG state measured at low temperature by 
25 37,39]. In our view, the bona fide PG state corre¬ 


sponds to the strict vanishing of both ft and S, the two 
parameters linked to the condensate. 

V. Energy phase diagram. — We have chosen a 


wide variety of tunneling DOS from the literature 18,28 
29p7p9pl]|55[|59 to study the energy scales, their doping 


dependence and the possible role of the pair-pair interac¬ 
tions. Five of these spectra labeled A-E in Fig. [5| have 
been selected for their regularity and the apparent stabil¬ 
ity of their parameters. In the upper panel, the QP DOS, 


using eq. (28), corresponds to underdoped, with larger 
A p ~ 50 meV, to overdoped, with smallest A p ~ 27 meV. 

In the lower panel, Fig. [5j we show the gap functions 
A k(Ek) used to precisely fit each spectrum of the upper 
panel. The experimental data are only shown for spectra 
A and E for clarity and, as in Fig. [4j only the SC part of 
the spectrum is considered. The quasiparticle peaks in the 
DOS are fairly sharp, reaching ^ 1.5 — 2 above the back¬ 
ground level. These are determined by the negative slope 
of the gap function <9A k{Ek)/dEk at Ek = A p . Higher 
QP peaks have sometimes been measured, in particular 
Ref. [59 , in good agreement with our model for large ft 
and smaller ao- 

Fig-i reveals the regularity of the dip position, given 
by the positive derivative of the gap-function near ^ 2 A p , 
which increases from the overdoped to underdoped sides. 
22|41|42 


As in Refs. 


these characteristics are found with¬ 


out a dynamical collective mode. It remains to be investi¬ 
gated whether similar features seen in the pnictide LiFeAs 
j|44 have the same origin. 

The parameters of these five spectra are then plotted 
in the T — p phase diagram, Fig. [6] The self-consistent 
gap A p follows a near-perfect linear trend decreasing 
from underdoped to overdoped sides and extrapolating 
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Fig. 6: Energy - doping (p) plot summarizing the dependence 
of the parameters in the present model. Points A-E correspond 
to the spectra of Fig. [5j The self-consistent gap A p is linear; 
Ao, <Jo and Edi P are decreasing with p but are slightly convex. 
At optimal doping e c , the pair-pair interaction f3 and the mini¬ 
gap 5 are maximal. Both (3 and ~ 7 5 follow the dome-shaped 
curve: 1.8 ksTc. 


Fig. 5: Upper panel: QP DOS using eq. ([28]) to fit a series of 
5 spectra A-E from various authors (A: [28] , B : 39], C : [55], 
D : [28], E: 129[) ordered according to gap width: A (over¬ 
doped), C (optimal), E (underdoped). (Data only shown in A 
and E.) Lower panel: Plot of the corresponding gap-functions 
A k(Ek) A-E with fitted parameters: gap width A p , minimum 
amplitude A 0 , and distribution width, ao- These parameters 
are plotted in the T — p diagram, Fig. [6] 


to zero at Q. To serve as a guide, the dome-shaped plot 
1.8 &bT c is shown. To a first approximation, for the dop¬ 
ing near the top of the T c -dome, the parameters Ao, 00 , 
and Edip are remarkably continuous and almost linear. 
This supports the conclusion that, excepting the highly- 
underdoped samples, the parameters deduced from the 
spectra of Fig. [5] are reliable. 


A closer look reveals the convex shape of Aq, with the 
consequence on e c = Ao — A p , as well as E^ip , also being 

evolves from below 2 A„ 


22 


convex. The dip position 18 
to above 2 A p for increasing p, as resolved by many tun¬ 
neling spectra. On the right-hand side, all curves seem to 
converge to a single point Q. 

Of major significance to the present model, the pair- 
pair interaction follows T c along with the attractive mini¬ 
gap 8, as plotted in the lower part of Fig. [6] We find /3 ~ 
1 .SksTc and 8 ~ .26&bT c with f3/8 ~ 7. Admittedly, 
/3 reveals some scatter, sometimes reaching ~ 2.2 ^T c . 
We note the loss of precision in the fits on the far under¬ 


doped side, where the coherence peaks are attenuated. In 
this regime where /? decreases sharply, the number of ex¬ 
cited pairs, related to cr, increases markedly, so that a 
practically joins cr 0 . The available data indicates that, at 
this extreme left-hand side of the p-diagram, both (3 and 
8 are vanishing. 

In our view, the overall scenario is as follows. A pair- 
pair interaction, governed by the coupling fyj allows for 
incoherent pairs, with distribution Po(A 2 ) to begin align¬ 
ing their energies, ultimately to the amplitude A p . Since 
the interaction is mediated by the quasiparticles, it pro¬ 
vides the mechanism for establishing long-range order. 
The final pair-pair interaction in the SC state, being pro¬ 
portional to the condensate density, /3(T) m /3 q N oc (T ), 
supports this hypothesis. At low-doping, the coherent 
state cannot be formed if the distance between pre-formed 
pairs is too large to be efficiently coupled, since the spa¬ 
tial inhomogeneity that breaks the coherence will result 
in a cut-off in the interaction. With increased doping, the 
interaction favors long-range order, the system becomes 
more homogeneous, and T c increases. At the same time, 
the pair-pair interaction weakens the pair self-energy A p 
(see fig.[6|, which decreases uniformly from under-doped 
to over-doped sides. Thus the SC state becomes weakened 
again. 

In the context of pre-formed pairs, the question of a 
non-BCS condensation mechanism arises. Indeed, in a 2D 
system, a conventional Bose condensation is not possible 
without an inter-particle interaction. We proposed that a 
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mini-gap S in the pair excitation spectrum allows for the 
condensation. The calculated condensate density N oc (T ) 
using this model follows very nicely the experimental data 
obtained both by single particle tunneling and by Joseph- 
son effect (fig. [ 3 ]). The temperature dependence thus has 
the properties of an order parameter; at the critical tem¬ 
perature, N oc (T c ) = 0 and we get back the incoherent 
pairs of the PG state. The mini-gap S is the direct mea¬ 
sure of the stability of the coherent state and, just like /?, 
it follows the T c versus p curve. 

Since S is the energy to excite a single pair with respect 
to the condensate, it can also be viewed as the energy 
needed to create a local defect in the phase. Such a defect 
might be the rotation of a local pair wave function with 
respect to the d-wave SC condensate. Without the mini¬ 
gap, there is no condensation at all and no coherent pair- 
pair interaction either (/3 = 0) - the two being intimately 
linked. 

In contrast to 6, which is the excitation energy of one 
pair, £ c = A 0 — A p represents the adiabatic condensation 
energy of all pairs from a PG state at zero-temperature 
to the SC state. This energy is of the same order of mag¬ 
nitude as p to which it is proportional, see eq. ( |Tt| ) , and 
therefore it too follows T c . The quantity e c also reveals 
the unconventional shift in pair binding-energy : from Aq, 
the average value in the (T = 0) PG state, to the smaller 
value A p in the SC state. This suggests that the pairs give 
up potential energy when the system acquires long-range 
order via the quasiparticles. This mechanism, important 
at low doping, becomes increasingly inefficient as doping 
becomes large and £ c is correspondingly small. Indeed, in 
the over-doped regime the behavior is dominated by the 
weakening of the superconductivity associated with the 
decrease of A p with p. Ultimately, all parameters vanish 
at a unique critical point Q. 

The optimal doping appears at the best compromise 
between the depairing associated with the decrease of A p , 
and the increase of the homogeneity and pair concentra¬ 
tion that stabilise the coherent state - both effects being 
driven by the pair-pair interaction. 

VI. Conclusion. — In this work we have shown that 
the phase transition to superconductivity from a pseudo¬ 
gap state of preformed pairs can be understood in rela¬ 
tively simple terms and few parameters. It consists of a 
Bose-like condensation with a small mini-gap S above A p 
which represents the minimum energy to remove a pair 
from the condensate. 

The low-level pair excitations above A p are thus critical 
for the condensation mechanism. No thermal QP excita¬ 
tions, a la BCS, are needed; only the pair excitation dis¬ 
tribution Po(Afc) is involved. With this mechanism, the 
properties of the condensate occupation number N oc (T ) 
are in good agreement with available experiments. More¬ 
over, the Fermi-level spectral gap A p remains constant in 
this temperature range. 

The problem is rendered complex by two additional 


properties of the pair-pair interaction • rsj 2/3(T)P 0 (A|). 
First, after examining a wide range of experimental spec¬ 
tra, we proposed that the excited states are coupled to 
quasiparticles A\ ^ E l k , which modifies the final gap 
equation, eq. (25). Once in the superconducting state, 
these QP excitations become well defined, revealing the 
characteristic peak-dip structures in the spectral function. 
The non-retarded energy-gap function described herein is 
thus the direct consequence of this QP-mediated interac¬ 
tion. 

Secondly, the pair-pair interaction also depends directly 
on the condensate: /3(T) = /3o N oc (T), which we identify 
as the order parameter. Not only does it vanish at T c , but 
it is directly linked to long-range order. This interaction 
varies spatially in the inhomogeneous superconductor, but 
also as a function of concentration p. In the latter case, 
the interaction is larger with increasing p which leads to 
a decrease in the self-consistent gap A p and consequently 
to the lowering of T c as well. Paradoxically, the very po¬ 
tential that is responsible for superconducting order even¬ 
tually, in the high-density limit, provokes its demise. 


REFERENCES 

[1] H.E. Stanley, Introduction to Phase Transitions and Criti¬ 
cal Phenomena (Oxford University Press, 1971) 

[2] S. Nazarenko, M. Onorato, and D. Proment, Phys. Rev. A 
90, 13624 (2014) 

[3] J. Bardeen, L. Cooper, J. Schrieffer, Phys. Rev. 108, 1175 
(1957). M. Tinkham, Introduction to Superconductivity, 
2nd Edition (McGraw-Hill, New York, 1996). 

[4] P. W. Anderson, Science 235, 1196 (1987) and refs, therein. 

[5] B. Edeggera, V. N. Muthukumar and C. Gros, Adv. Phys. 

56, 927 (2007). 

[6] J. Maly, B. Janko, and K. Levin, Phys. Rev. B 59, 1354, 
(1999) ; V. Loktev, R. Quick, S. Shaparov, Phys. Rep. 349, 
1 (2001), and refs, therein. 

[7] N.-C. Yeh and A. D. Beyer, Int. J. Mod. Phys. B 23, 4543 
(2009) and refs, therein. 

[81 N. Andrenacci, M. Capezzali and H. Beck, Eur. Phys. J. B 
53, 417 (2006) 

[9] V. Emery and S. Kivelson, Nature 374, 434 (1995) 

[10] Spectral function M. R. Norman, M. Randeria, B. Janko, 
J. C. Campuzano, Phys. Rev. B 61, 14742 (2000). 

[11] A. Chubukov and D. Morr, Phys. Rev. Lett. 81, 4716 
(1998). 

[12] M. Eschrig, Adv. Phys. 55, 47 (2006). 

[13] M. Norman, M. Randeria, H. Ding et al., Phys. Rev. B 

57, R11093 (1998); J.C. Campuzano et al., Phys. Rev. Lett. 
83, 3709 (1999). 

[14] J.C. Campuzano et al., Phys. Rev. Lett. 83, 3709 (1999). 

[15] M. R. Norman, H. Ding, Phys. Rev. B 57, R11089 (1998). 

[16] N. Jenkins, Y. Fasanol, C. Berthod, I. Maggio-Aprile, A. 
Piriou, E. Giannini, B. W. Hoogenboom, C. Hess, T. Cren, 
and 0. Fischer, Phys. Rev. Lett. 103, 227001 (2009). 

[17] G. Levy de Castro, Ch. Berthod, A. Piriou, E. Giannini, 
and 0. Fischer, Phys. Rev. Lett. 101, 267004 (2008) 

[18] J. F. Zasadzinski, L. Ozyuzer, N. Miyakawa, K. E. Gray, 
D. G. Hinks, and C. Kendziora, Phys. Rev. Lett. 87, 67005 


P-9 




W. Sacks, A. Mauger & Y. Noat 


( 2001 ). 

[19] W.Y. Liang, J. Phys. Condens. Matter 10, 11365 (1998); 
T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999) 

[20] M. Le Tacon, A. Sacuto et al., Nature Physics 2, 537 
(2006). 

[21] C. Huscroft and R. T. Scalettar, Phys. Rev. Lett. 81, 2775 
(1998). 

[22] W. Sacks, T. Cren, D. Roditchev and B. Dougot, Phys. 
Rev. B 74, 174517 (2006). 

[23] J. Tallon and J. Loram, Physica C 349, 53 (2001). 

[24] S. Pan, J. P. O’Neal, R. L. Badzey et ah, Nature 413, 282 

( 2001 ). 

[25] T. Cren, D. Roditchev, W. Sacks, J. Klein, J.-B. Moussy, 
C. Deville-Cavellin, and M.Lagues, Phys. Rev. Lett. 84, 
147 (2000) 

[26] T. Cren, D. Roditchev, W. Sacks et ah, Europhys. Lett., 
54 (1), 84 (2001). 

[27] C. Howald, P. Fournier and A. Kapitulnik, Phys. Rev. B 
64, 100504 (2001). 

[28] K. McElroy, Jinho Lee, J. Slezak, D.-H. Lee, H. Eisaki, S. 
Uchida, J.C. Davis, Science 309, 1048 (2005) 

[29] A. Sugimoto, S. Kashiwaya, H. Eisaki, H. Kashiwaya, H. 
Tsuchiura, Y. Tanaka, K. Fujita, and S. Uchida, Phys. Rev. 
B 74, 94503 (2006) 

[30] J. Hoffmann, E. Hudson, K. Lang, H. Eisaki, B. Uchida 
and J.C. Davis, Physica C 388, 703 (2003) 

[31] I. Zeljkovic, E. Main, et ah, Nature Materials 11 585 

( 2012 ). 

[32] H. Alloul, T. Ohno and P. Mendels, Phys. Rev. Lett. 100, 
226802 (1989). 

[33] H. Ding, J.C. Campuzano, M. Norman et ah, J. Phys. 
Chem. Solids 59, 1888 (1998). 

[34] M. Franz and A.J. Millis, Phys. Rev. B 58, 14572 (1998). 

( 2000 ). 

[35] A. Yazdani et ah, J. Phys.: Condens. Matter 21, 164214 
(2009). 

[36] A. Matsuda, S. Sugita, T. Fujii, T. Watanabe, J. Phys. 
Chem. Sol. 62, 65 (2001) 

[37] Ch. Renner et ah, Phys. Rev. Lett. 80, 149 (1998). 

[38] Ch. Renner et ah, Phys. Rev. Lett. 80, 3606 (1998). 

[39] S. H. Pan, E. W. Hudson, A. K. Gupta, K.-W. Ng, H. 
Eisaki, S. Uchida, and J. C. Davis, Phys. Rev. Lett. 85, 
1536 (2000). 

[40] 0. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod and 
C. Renner, Rev. Mod. Phys. 79, 353 (2007). 

[41] T. Cren, D. Roditchev, W. Sacks and J. Klein, Europhys. 
Lett. 52 (1), 203 (2000). 

[42] W. Sacks, A. Mauger, Y. Noat, Physica C 503, 14 (2014) 

[43] M. Norman et ah, Phys. Rev. Lett. 79, 3506 (1997); M. 
Eschrig and M. Norman, Phys. Rev. Lett. 85, 3261 (2000). 

[44] Shun Chi, S. Grothe, Ruixing Liang, P. Dosanjh, W. N. 
Hardy, S. A. Burke, D. A. Bonn, and Y. Pennec, Phys. 
Rev. Lett. 109, 087002 (2012). 

[45] O. Ahmadi, L. Coffey, and J. F. Zasadzinski, N. Miyakawa, 
L. Ozyuzer, Phys. Rev. Lett. 106, 167005 (2011) ; and 
references therein. 

[46] P. Romano, L. Ozyuzer, Z. Yusof, C. Kurter, and J.F. 
Zazadzinski, Phys. Rev. B 73, 92514 (2006). 

[47] A. Ghosal, M. Randeria, and N. Trivedi, Phys. Rev. B 65, 
014501 (2001). 

[48] W. A. Atkinson, Phys. Rev. B 71, 24516 (2005). 

[49] W. A. Atkinson, P. J. Hirschfeld, and Lingyin Zhu, Phys. 


Rev. B 68, 054501 (2003). 

[50] B. M. Andersen, S. Graser, and P. J. Hirschfeld, Phys. 
Rev. Lett. 105, 147002 (2010). 

[51] D. Valdez-Balderas & D. Stroud, Phys. Rev. B 74, 174507 
(2006) 

[52] M. Tinkham, Introduction to Superconductivity, 2nd Edi¬ 
tion (McGraw-Hill, New York, 1996). 

[53] R.C. Dynes, V. Narayanamurti, J.P. Garno, Phys. Rev. 
Lett. 41 (1978) 1509. 

[54] A. Fetter and J. Walecka, Quantum theory of many parti¬ 
cle systems , (McGraw-Hill, New York, 1971). 

[55] N. Miyakawa, J. Zasadzinski, L. Ozyuzer, P. Guptasarma, 
D. Hinks, C. Kendziora, and K.E. Gray, Phys. Rev. Lett. 
83, 1018 (2001). 

[56] J.R. Schrieffer, Rev. Mod. Phys. 200, (1964), and 
J.R.Schrieffer, Theory of Superconductivity , (W.A. Ben¬ 
jamin, New York, 1964). 

[57] P. Mallet, D. Roditchev, W. Sacks, D. Defourneau, and J. 
Klein, Phys. Rev. B 54, 13324 (1996). 

[58] W. Hoogenboom, C. Berthod, M. Peter, 0. Fischer, and 
A. A. Kordyuk, Phys. Rev. B 67, 224502 (2003). 

[59] A. Fang et ah, Phys. Rev. Lett. 96, 017007 (2006). 


p-10 



